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ABSTRACT 


We  evaluate  the  bit  error  probability  of  coded  lightwave  systems  employing 
direct  detection  with  00K  modulation,  and  coherent  optical  systems  with  FSK 
modulation  and  noncoherent  detection.  Block  codes,  convolutional  codes,  and  con¬ 
catenated  codes  are  investigated.  For  direct  detection  systems,  both  hard  and  soft 
decoding  are  considered.  Only  hard  decoding  is  employed  in  coherent  optical  sys¬ 
tems. 
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I.  INTRODUCTION 


The  benefits  of  forward  error-correction  coding  for  radio  communications  in  an 
additive  white  Gaussian  noise  (AWGN)  channel  are  well  documented  [1-3].  With 
the  state-of-the-art  technology  of  Application  Specific  Integrated  Circuits  (ASIC), 
both  block  coding  and  convolutional  coding  can  be  implemented  at  very  high  bit 
rates.  Block  coding  can  be  implemented  at  Gb/s  data  rates  and  convolutional  coding 
at  hundreds  of  Mb/s  data  rates.  This  opens  the  possibility  for  employing  FEC  in 
lightwave  systems  to  compensate  for  the  loss  in  the  link  budget  or  to  relax  the  laser 
linewidth  requirement  in  coherent  optical  systems  [4-5]. 

In  this  study,  we  explore  the  benefit  of  FEC  for  direct  detection  lightwave 
systems  employing  on-off-keying  (00K)  modulation,  and  coherent  optical  systems 
employing  frequency-shift  keying  (FSK)  with  envelope  detection  corrupted  by  phase 
noise.  Our  work  was  derived  explicitly  from  [19]. 

Of  practical  importance  is  the  use  of  concatenated  codes.  Concatenated  codes 
allow  the  use  of  two  relatively  short  codes  with  modest  error-correcting  capability, 
one  as  the  inner  code  and  the  other  as  the  outer  code  to  achieve  a  large  coding  gain. 
Short  codes  with  modest  error-correcting  capability  simplify  decoder  design  at  very 
high  data  rates  and  thus  reduce  the  implementation  cost.  Concatenated  codes  also 
allow  the  use  of  a  short  constraint  length  convolutional  code  as  the  inner  code  with 
soft  decoding  and  a  short  Reed-Solomon  (RS)  code  with  modest  error-correcting 
capability  as  the  outer  code.  This  scheme  can  achieve  a  large  coding  gain  which 
otherwise  can  only  be  achieved  with  a  long  constraint  length  convolutional  code  or 
a  long  RS  code  with  large  error- correcting  capability. 
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II.  BIT  ERROR  PROBABILITY  OF  DIRECT 
DETECTION  LIGHTWAVE  SYSTEMS 


A.  OOK  MODULATION 

The  direct  detection  OOK  receiver  model  to  be  analyzed  is  given  in  Fig.  1.  The 
photodiode  detector  has  a  responsivity  %  =  qq/hj  (A/W)  where  q  is  the  quantum 
efficiency  which  is  near  unity  for  p  —  i  —  n  diodes,  q  =  1.6  x  10-19  C  is  the  electron 
charge,  h  =  6.626  x  lO-34  J  S  is  Planck’s  constant,  and  /  is  the  laser  frequency. 
The  current  produced  by  the  photodiode  is  amplified  by  a  low  noise  amplifier  and 
integrated  by  the  integrator.  With  soft  decoding,  the  analog  samples  from  the 
integrator  are  processed  by  the  soft  decoder.  With  hard  decoding,  the  integrator's 
samples  are  first  detected  by  a  threshold  detector  to  provide  hard  quantized  samples 
to  the  hard  decoder. 

In  practice,  the  shot  noise  and  dark  currents  produced  by  the  photodiode 
are  negligible  compared  to  the  thermal  noise  current  generated  by  the  low  noise 
amplifier.  Therefore,  they  are  ignored  in  our  analysis.  The  complex  envelope  of  the 
OOK  lightwave  signal  incident  upon  the  photodiode  detection  during  a  coded  bit 
interval  T  is 

a,-(0  =  y/p  b:d9W  i  =  0, 1  ;  0  <  t  <  T  (1) 


OOK 


•DATA 


Noise 

n(<) 


Figure  1:  Coded  OOK  direct  detection  receiver  structure. 
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where  P  is  peak  power,  Oq  —  \Jpj{  1  -f  p),  b\  =  ^1/(1  +  p)  and  p  is  the  laser  extinc¬ 
tion  ratio.  The  output  of  the  integrator  is 

Yt  =  K  F  \si(t)\2dt+  f  n(t)dt  (2) 

Jo  Jo 

where  n(t)  is  the  thermal  noise  current  with  spectral  density  No  (A2/Hz).  Substi¬ 
tuting  (1)  into  (2)  and  denote 

N  =  [T  n(t)dt  (3) 

Jo 


we  obtain 


Yt  =  HPTb]  +  N 


(4) 


Note  that  N  is  a  Gaussian  random  variable  with  zero  mean  and  variance  cr2  =  NqT. 
We  assume  that  the  soft  decoding  is  done  by  a  maximum  likelihood  soft  decoder 
using  the  analog  (unquantized)  samples  Yi  in  (4).  For  hard  decoding,  a  bit  decision 
at  the  output  of  the  threshold  detector  employing  a  threshold  setting  a  is  produced 
with  channel  transition  probability  as  follows: 


P  =  iPr{yi>a}  +  ij>r{y,<a} 


(5) 


where  Q(x)  is  give  by 

(6) 

The  minimum  value  of  p  at  the  optimum  threshold  setting  aopt  =  72.FT’(62  +  &o)/2 
is 


P  + 


„  (HPT®  ~  >S)\ 
2V7W  ) 


=  q(~SNR ,/f) 


(7) 
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where 


SNR  =  TlP{b\  -  6o)y  (8) 

is  the  signal- to- noise  ratio  at  the  threshold  detector  input,  r  is  the  code  rate,  and 
Ti  is  the  bit  interval. 

B.  CODING 

Assuming  that  the  decoding  is  accomplished  by  a  maximum  likelihood  decoder, 
then  the  bit  error  probability  of  the  decoded  00K  signals  employing  block  code 
(n,  k)  is  upper  bounded  by 

a  <  E  o) 

where  rij  is  the  number  of  code  words  of  weight  j,  and  P}  is  the  probability  that  the 
distance  between  the  received  sequence  and  a  weight- j  code  word  is  less  than  the 
distance  to  an  all-zero  code  word 

For  a  convolutional  code  with  rate  k/n,  the  bit  error  probability  is  upper 
bounded  by  [1-3] 

a  <  j  £  *>p>  o») 

where  u>j  is  the  total  information  weight  of  all  code  paths  of  weight  j,  dj  is  the  free 
distance  of  the  code,  and  Pi  is  the  probability  that  the  all- zero  path  is  eliminated 
by  a  path  of  weight  j  merging  with  it. 

With  hard  decoding,  the  probability  P;  in  (9)  and  (10)  is  given  by 

5Z  ( { jpii-pr.  i°dd 

,=(j+l)/2  \  / 

P,  =  j  -  (11) 

\  (  // 2  ) 7/72(1  “  p)j/2  +  ^  (  \  )p'(1  ~  3  even 
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where  p  is  the  channel  transition  probability  in  (7).  Alternatively  P0  in  (9)  can  be 
calculated  by  the  following  upper  bound 


a<  E  —  <12> 

j=<+ 1  n  \  3  ) 

where  t  is  the  error  correcting  capability  of  the  block  code. 

For  soft  decoding,  the  decoder  selects  the  code  word  (or  code  path)  at  minimum 
Euclidean  distance  from  the  received  sequence.  Thus,  Pj  in  (9)— ( 10)  is  the  probabil¬ 
ity  that  the  Euclidean  distance  between  the  received  sequence  Y  and  weigh  t-j  code 
word  Xj  is  less  than  the  Euclidean  distance  to  the  assumed  transmitted  all-zero 
code  word  A'0.  Thus 

P,  =  Pr  {||V  -  A',111  <  ||V  -  A'oll1}  (13) 

where  Xj  is  a  vector  with  j  components,  each  of  value  HPTb\,  represented  bit  ones, 
and  the  remaining  components,  each  of  value  KPTb\ ,  represented  bit  zeros;  the  vec¬ 
tor  Xq  represents  the  assumed  transmitted  all-zero  sequence  with  each  component 
assuming  the  value  TtPTb q.  The  received  vector  is  Y  —  Xo+  N  where  N  is  the  zero 
mean  Gaussian  noise  vector,  each  component  has  variance  tr2  =  N0T. 

We  note  that  by  using  Y  =  X0  +  N  in  (12)  we  have 

Pi  =  Pr{||Ar-(A,-A„)||J<||)V||!} 

=  Pr  {l|(^llJ  +  IIA,  -  Aof  -  2 Re  {(N, X,  -  A„)}  <  ||JVf } 

=  Pr{2Re  {(N,X,  -  A„))  >  ||A,  -  A„||3} 

-  (H) 

where  (N,Xj  —  X0)  is  the  inner  product  of  the  two  vectors  N  and  Xj  —  Xo ■  We 
remark  that  Re  {(iV,  (Xj  -  A0)/||X;  -  A'o||}}  is  the  sum  of  zero  mean  Gaussian 


random  variables,  each  with  variance  a 2  =  N0T,  and  (Xj  —  Xo)l\\Xj  —  Xo||  is  a  unit 
length  vector.  Thus,  it  is  a  zero  mean  Gaussian  variable  with  variance  a2  =  N0T. 
Hence 


_  n(KPn%-%)sri\ 

~  yl  2  VAST  ) 

= 

=  q(\snr^) 


(15) 


where  r  is  the  code  rate  and  SNR  is  given  in  (8). 

The  above  results  are  applied  to  various  codes  to  study  their  effectiveness. 
Figure  2  shows  the  bit  error  probability  Pb  versus  signal-to-noise  ratio  SNR  for 
the  (7,  4)  Hamming  code  and  the  concatenated  (7,  4)/(15,  11)  Hamming/Reed- 
Solomon  codes  with  the  Hamming  code  as  the  inner  code.  With  hard  decoding  the 
Hamming  code  provides  very  small  coding  gain.  On  the  other  hand,  the  coding  gain 
at  Pb  =  10~15  is  about  2.3  dB  for  the  concatenated  codes  with  soft  decoding  for  the 
Hamming  code.  The  coding  gain  at  Pb  =  10-15  is  almost  3  dB  as  shown  in  Fig.  3 
for  the  concatenated  (31,  27)/(27,  23)  RS/RS  with  hard  decoding  where  the  (27, 
23)  is  the  shortened  outer  code. 

Figure  4  shows  the  performance  of  a  convolutional  code  with  rate  3/4  and 
constraint  length  v  =  2  and  v  =  5.  The  shorter  constraint  length  v  =  2  is  more 
suitable  for  high  data  rates.  The  coding  gain  at  Pb  =  10~15  is  about  1.5-dB  for  v  =  2 
and  2.5  dB  for  v  =  5.  Large  coding  gain  can  be  obtained  with  the  concatenated  rate 
3/4  convolutional/(31,  27)  RS  codes  with  the  convolutional  code  as  the  inner  code 
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as  in  Fig.  5.  A  coding  gain  of  about  2.8  dB  and  4  dB  can  be  achieved  at  Pb  =  10-15 
with  soft  decoding  for  the  convolutional  code  with  v  =  2  and  v  =  5,  respectively. 


S 
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Figure  3:  Coded  Pt,  versus  SNR  for  direct  detection  systems  employing  (31, 
27 ) / (27,  23)  Reed-Solomon/Reed-Solomon  codes. 
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III.  BIT  ERROR  PROBABILITY  OF  COHERENT 
OPTICAL  FSK  SYSTEMS  WITH  NONCOHERENT 

DETECTION 


Heterodyne  and  homodyne  (coherent)  lightwave  systems  are  largely  affected 
by  the  combined  laser  phase  noise  of  the  source  and  local  laser  oscillators.  The 
performance  deteriorates  as  the  linewidth  bit  time  0Tt,  increases  where  0  (Hz)  is 
the  combined  laser  linewidth  of  the  source  and  local  lasers  and  T*.  (s)  is  the  bit 
time  [6-10].  Coherently  detected  signals  suffer  much  more  from  the  phase  noise 
than  noncoherently  detected  signals.  The  phase  noise  which  can  be  modeled  as  a 
Brownian  motion  process  causes  signal  fading  and  broadens  the  signal  spectrum. 
Both  effects  can  be  compensated  to  a  large  extent  by  using  a  wider  predetection 
filter  to  pass  more  of  the  signal  power  and  an  averaging  postdetection  filter  to 
provide  time  diversity  and  thus  alleviating  the  signal  fading  effect.  Such  a  technique 
was  investigated  in  [11-12]  for  heterodyne  lightwave  systems  employing  ASK  and 
wideband  FSK  with  noncoherent  detection,  and  in  [13]  for  homodyne  lightwave 
systems.  The  use  of  time  diversity  is  equivalent  to  the  use  of  simple  repetition 
codes.  Since  simple  repetition  codes  can  provide  a  good  improvement  as  reported  in 
[11-12],  one  would  naturally  expect  an  even  larger  improvement  when  more  powerful 
and  efficient  codes  are  used. 

In  this  part  we  investigate  the  performance  of  heterodyne  lightwave  systems 
using  wideband  orthogonal  FSK  with  noncoherent  detection  and  error- correction 
codes.  We  will  show  that  coding  can  provide  a  substantial  improvement  in  the  face 
of  the  phase  noise  and  shot  noise  modeled  as  a  zero-mean  Gaussian  process  with 
power  spectral  density  (PSD)  ./Vo/2  (W/Hz).  Furthermore,  coding  can  relax  the  0T^ 


requirement  a  great  deal  thus  enabling  the  use  of  inexpensive  lasers  with  a  much 

larger  linewidth  than  that  required  by  an  uncoded  system.  For  total  generality  we 

consider  both  time  diversity  and  coding  for  the  system  under  investigation  as  shown 

in  Fig.  6.  We  assume  that  the  frequency  spacing  j/0  —  fi\  is  much  larger  than  1/T' 

where  V  is  the  integration  time  and  is  also  the  sampling  time  at  the  summer  outputs. 

The  time  T'  is  selected  such  that  the  code  symbol  time  T  =  rT*  =  MT'  where  M  is 

a  positive  integer  and  r  is  the  code  rate.  Time  diversity  of  the  postdetection  filter  is 

modeled  by  summing  M  independent  samples  of  a  coded  bit  after  envelope  detection. 

The  output  samples  are  fed  to  the  decoder  which  can  make  a  hard  decision  or  soft 

decision  on  the  received  codewords.  Equivalently,  each  branch  of  the  demodulator 

in  Fig.  6  can  be  replaced  by  a  bandpass  integrator  followed  by  a  square-law  envelope 

detector  [12].  In  the  following  section  we  derive  the  channel  transition  probability. 

We  provide  a  simple  derivation  as  an  alternative  to  the  one  in  [12]. 

For  mathematical  convenience  we  adopt  the  envelope  notation  of  a  real  signal. 

Thus,  the  summer  outputs  of  the  signal-present  upper  branch  and  the  signal-absent 

lower  branch  in  Fig.  6  with  integration  time  T'  are 

M 

y<)  =  53 1-^*  Ucok  +  jn»o*i2  (i6) 

k- 1 

M 

Y\  -  X)  K»*  +  >n»i*l2  (17) 

*=i 

where  Zk  is  given  by 

Zk  =  ~  ["  e’mdt  (18) 

and  all  the  in-phase  and  quadrature  noise  samples  ncik ,  nsik(i  =  0, 1;  k  =  1,2,-  *  - ,  M) 
are  independent  and  identically  distributed  (i.i.d.)  zero-mean  Gaussian  random 
variables  (R.V.)  with  variance  NoT/4.  The  laser  phase  noise  process  0(f)  can  be 
characterized  as  a  zero-mean  Gaussian  process  with  variance  2-nfit  [16],  where  /?  is 
the  combined  linewidth  of  the  source  and  local  laser  oscillators. 
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Figure  6:  Coded  coherent  optical  FSK  receiver  structure. 

Since  Vj  is  the  sum  of  squares  of  2 M  i.i.d.  Gaussian  R.V.s  with  zero  mean  and 
variance  a 2  —  N0T'/4,  Y\  is  chi-square  distributed  with  the  following  probability 
density  function  (pdf)  [14] 

^  ~  (M  -  l)!2Mo,2M^  l£  n/2°  ’  yi  “  ° 

On  the  other  hand,  Vo  is  the  sum  of  squares  of  M  independent  Gaussian  R.V.s  with 
variance  a 2  and  mean  and  M  independent  Gaussian  R.V.s  with  the  same 

variance  <r2  but  with  mean  /m{JYfe}.  Here  Re{,Yt}  and  Im{Xk}  denote  the  real  and 
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imaginary  parts  of  A*,  respectively.  Thus,  Yo  is  non-central  chi-square  distributed 
with  the  following  conditional  pdf  [14] 

Vo>0  (20) 

where  lsi~ i(*)  is  the  ( M  —  l)th-order  modified  Bessel  function  of  the  first  kind,  and 
v  is  the  value  assumed  by  the  random  variable  V  defined  as 

V=£|Z*|2  (21) 

fcssl 

The  conditional  channel  transition  probability  of  error  is  given  by 


p(v)  =  Pr  >  Volf} 

noo 

/v,(yi)/yo(yo!w)dyidy0 

0 


The  inner  integral  can  be  evaluated  as  follows: 


Substituting  (20)  and  (23)  into  (21),  we  obtain 

pH  =  %  h 

where  w  is  the  value  assumed  by  the  random  variable  W  as  defined  as 

V  i  M 

Using  the  identity  [15] 

xv^+neix Iv(2ay/x  )dx  = 


where 


K(*)  =  ±(-ir(  n+v)s^ 

ho  \n-m  )  m! 


(22) 

(23) 

(24) 

(24) 

(26) 

(27) 
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in  (24)  we  obtain 


p{w)  = 


e~w/2  *^llk{k  +  M-l\  {w/2)r 


2m 
s-u</2  M- 1 


M- 1  ,  *  / 

E  jrE  ( 

fc=0  c  m=0  \ 


k  —  m 


ml 


2m 


^  (w/2)n  ^  / 
n=0  ”•  i=n  V 


A:  +  M  —  1  \  1 


k  —  n 


2k 


Define 


M 


SNR  =  A2Tb/2N0  ,  y  =  '£)  Xk{rc/M) 


k=  1 


where  r  is  the  code  rate,  c  =  2irflTb  and  Xk  is  given  by 


Xk(rc/M)  = 


[k 

Jk-l 


(28) 


(29) 


with  rj;(t)  as  a  zero  mean  Gaussian  process  of  variance  t ;  we  can  express  the  condi¬ 
tional  channel  transition  probability  in  (28)  as  follows: 


p{y ) 


=  e~(SNRW2M  *=»  [( SNR)ry/2M]n  ^  f  k  +  M  -  l  \  1. 

fc=n  V 


2M 


n=0 


n\ 


k  —  n 


¥  <»> 


The  channel  transition  probability  is  given  by 


P=E{p(y)}=[  p{y)fY{y)dy  (31) 

where  fy(y)  is  the  probability  density  function  of  y  and  can  be  calculated  by  the 
method  given  in  [12].  Figure  7  shows  fy(y)  for  M  =  2,  c  =  1,  and  c  =  4  for  the 
uncoded  case. 

The  transition  probability  p  in  (31)  is  applied  to  (9)— (12)  to  evaluate  the 
effectiveness  of  various  codes  for  coherent  optical  systems.  Figure  8  shows  the  bit 
error  probability  Pb  versus  the  signal-to-noise  ratio  SNR  for  M  =  1  and-c  =  1  (small 
phase  noise)  with  hard  decoding  (15, 11)  and  (31,  27)  RS  codes  and  (31,  27)/(27,  23) 
concatenated  RS/RS  codes.  It  is  seen  that  coding  provides  an  excellent  performance 
as  compared  to  an  uncoded  system,  especially  for  the  concatenated  code.  When 
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M  =  2,  that  is,  the  bandwidth  of  the  integrators  is  twice  the  coded  rate,  more 
signal  power  passes  through  the  demodulator  and  the  performance  of  an  uncoded 
system  as  shown  in  Fig.  9  improves  remarkably  over  an  uncoded  system  for  M  —  1 
in  Fig.  8.  In  this  situation  the  improvement  in  performance  when  RS  codes  are 
employed  is  less  than  that  in  Fig.  8  but  the  improvement  is  still  considerable, 
especially  for  the  concatenated  (31,  27) / (27,  23)  RS/RS  codes.  When  the  phase 
noise  is  large  as  shown  in  Fig.  10  for  c  =  4,  the  use  of  both  RS  codes  and  larger 
integrator  bandwidth  M  =  2  achieves  considerable  improvement  for  the  (15,  11)  and 
(31,  27)  RS  codes  and  the  concatenated  (31,  27)/(27,  23)  RS/RS  codes.  Figures  11 
and  12  show  the  performance  of  rate  3/4  convolutional  codes  with  hard  decoding  for 
M  =  2,  and  c  =  1  and  c  =  4,  respectively.  It  is  observed  that  when  the  constraint 
length  is  small,  v  =  2,  the  improvement  is  small.  The  performance  of  rate  3/4. 
v  —  5  convolutional  codes  is  slightly  less  than  that  of  RS  codes  (15,  11)  or  (31,  27). 
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9 


Probability  of  Bit  Error 


Figure  8:  Coded  P \  versus  SNR  for  coherent  optical  systems  employing  the  (15, 
11)  and  (31,  27)  Reed-Solomon  codes  and  the  concatenated  (31,  27)/ (27,  23)  Reed- 
Solomon/ Reed- Solomon  codes  with  M  =  1  and  c—  1. 


20 


Probability  of  Bit  Error 


SNR(dB) 


Figure  9:  Coded  Pj,  versus  SNR  for  coherent  optical  systems  employing  the  (15, 
111  and  (31,  27)  Reed-Solomon  codes  and  the  concatenated  (31,  27)/(27,  23)  Reed- 
Solomon/ Reed- Solomon  codes  with  M  =  2  and  c  =  1. 
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Figure  10:  Coded  Pt,  versus  SNR  for  coherent  optical  systems  employing  the  (15, 
11)  and  (31,  27)  Reed-Solomon  codes  and  the  concatenated  (31,  27)/(27,  23)  Reed- 
Solomon/Reed-Solomon  codes  with  M  =  2  and  c  =  4. 


SNR  (dB) 


Figure  11:  Coded  Pt,  versus  SNR  for  coherent  optical  systems  employing  hard 
decoding  convolutional  code  with  rate  3/4  and  constraint  length  2,  M  =  2,  c  =  1 
and  c  =  4. 
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IV.  CONCLUSIONS 


From  our  investigaion  we  conclude  that  error  correction  coding  can  improve 
the  performance  of  both  direct  detection  and  coherent  optical  systems.  The  im¬ 
provement  is  considerable  when  concatenated  codes  are  used.  From  the  practical 
point  of  view,  short  RS  codes  with  modest  error  correcting  capability  can  be  con¬ 
catenated  to  provide  large  coding  gain  and  at  the  same  time  can  be  implemented  at 
high  data  rates.  Convolutional  codes  with  soft  decision  also  provides  good  coding 
gain  at  moderate  data  rates  for  direct  detection  systems.  It  remains  to  be  seen 
how  much  coding  gain  can  be  obtained  with  soft  decision  convolutional  codes  for 
coherent  optical  systems. 

Our  contributions  in  this  thesis  are  the  numerical  results  obtained  from  the 
analytical  work  in  [19]  as  shown  in  equations  (1-31). 
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APPENDIX 


The  following  MATHEMATICA  (Stephen  Wolfram,  1993)  and  MATLAB  (The  Math 
Works,  1993)  codes  offer  the  formulations  for: 


The  convoh  aon  method  for  the  determination  of  the  probability  density  func¬ 
tion 


The  computation  of  bit  error  probability  for  direct  detection  and  coherent 
detection 

1.  Obtaining  of  pdf  of  y  by  M-fold  convolution 


xl*[ . ]  . Convolution  method  for  finding  pdf  of  y 

. Input  of  pdf  of  x 

x2*xl;  . . . xl,x2  are  i.i.d 

x3*0: 0.02:2; 

yl*0.02*conv(xl,x2) ;  . 0.02  is  integer  as  interval 

. 2nd  convolution 


2.  Bit  error  probability  for  direct  detection 
2.1  Uncoding  case 


function  y*q(x)  . Definition  of  q(x)  in  MATLAB 

y*0.5*erfc(x/sqrt(2))  . Use  error  function  defaulted  in 

MATLAB 

snrdb»6 .5:0.02:12.1; 
snr»10 . * (snrdb*0 . 1) ; 

pu»q(snr*0.5) ;  . Uncode  case 

2.2  Concatenated  convolutional  hard  and  soft/RS  coding 


Ml«input(' input  inner  convolution  code  order  as  n,k,v,df ,vd,vd+l . . . ') ; 

. Need  to  input 

two  matrix, one  for  RS  code(outter) , 
one  for  convolution  (inner) 

M2« input (’input  outter  RS  code  as  n,k’); 
nl«Ml(l);  kl«Ml(2);  n2«M2(l);  k2«M2(2); 

rl«kl/nl;  r2*k2/n2; 
m2«log2(n2+l) ;  t2«(n2-k2)*0.5; 
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r«rl*r2; 

p«q(snr*r“0.5*0.5) ;  . Equation  (7) 

pjl*[3;  . Solve  five  pj  values 

Hard  case 

pj2*[3 ;  .  Soft  case 

for  d*0:4 
pj»0; 

j“Ml (4)+d; 

if  rem(j,2)*«l  . Check  odd  or  even  case 

for  i*(j+l)/2:j 

pj«pj+bino(j ,i)*p.“i.*(l-p) .‘(j-i) ; 
end 


else 

for  i»j/2+l:j 

pj*pj+bino(j,i)*p.*i.*(l-p) .‘(j-i); 
end 

pj*pj+0.5*bino(j ,0.5*j)*p.“(0.5*j) .*(l-p) .  “(0.5*j) ; 

end 

pjl*[pjl;pj] ;  .  Hard  case's  pj 

PPjjm<l(  (r*j)“0.5*snr*0.5) ;  . Soft  case's  pj 

pj2«[pj2;ppjj3; 

end 

pb_conv=0;  . Solve  for  pb  for  hard  case 

ppbb_conv*0;  . For  soft  case 

for  h*l:5  . 5  is  from  vd  to  vd+4 

pb_  conv«pb_  conv+M 1 (h+4) *p  j 1 (h, : ) ; 

ppbb_conv«ppbb_conv+Ml (h+4) *pj  2 (h , : ) ; 

end 

pb_conv=pb_conv/kl ; 
ppbb  _  conv=ppbb_  conv/kl ; 

ps«l-(l-pb_conv) . “m2;  . Solve  for  ps 

ppss*l-(l-ppbb_conv) .“m2; 

pb*0;  ppbb*0;  . Solve  for  pb 

for  j*t2+l:n2 

pb*pb+(j+t2)/n2*bino(n2, j)*ps . “j .*(l-ps) . “(n2-j) ; 
ppbb»ppbb+(j+t2)/n2*bino(n2,j)*ppss.“j .*(l-ppss) .“(n2-j) ; 
end 

pb*pb*(n2+l)/(2*n2) ; 
ppbb»ppbb*(n2+l)/(2*n2) ; 

ss*linspace(l ,len,51) ;  . Plot  statement 

semilogy  (snrdb(ss) ,  pu(ss) , ,snrdb,pb, ' — ' ,snrdb,ppbb) 
axis([snrdb(l)  snrdb(len)  10“ (-15)  1]) 


3.  Bit  error  probability  for  coherent  detection 
3.1  Uncoding  case 


> 


yii-< 

M-2; 

S-55; 
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Input  pdf  value  of  y 
Wanted  SNR 


yl-Table[y,{y, 0.01761,1.9,0.02)] ; 
pl*Table [Sum  C ( ( (S*y ) / (2*M) ) “n) / (n ! ) 

♦Sum [Binomial [M+k-1 ,k-n3 
*(2“(-k)) ,{k,n,M-l)j  ,{n,0,M-l)3 
♦Exp [-(S*y/ (2*M))3/ (2“M) , 

{y, 0.01761, 1.9, 0.02)3; 

p2*Table [{y 1 [ [i] 3 ,pl [ [i] ] *yll [ [i3  3  > , 

{i,l, Length [yl3)3 ; 

p3*Interpolation[p2,  InterpolationOrder  ->  13 

.  Connect  each  point 

p4*NIntegrate[p3[y3 ,{y, 0.01761, 1 .89761)3 


3.2  concatenated  RS/RS  coding 


yll«{ . };  . For  coding  case’s  pdf  of  y  value 

For[M*l;  g*l;  nl»31;  n2*27;  kl*27;  k2=23;  rcl«kl/nl; 

rc2*k2/n2;  Sl*10;  . Inputs  of  SNR  and  code  rate 

S«Sl*rcl*rc2,  S  +«5, 
y l*Table [y , {y , 0 . 04283 , 0 . 97772 , 0 . 02)3 ; 
pl*Table [Sum [ ( ( (S*y) / (2*M) ) “n) / (n ! ) 

♦Sum [Binomial [M+k-1 ,k-n3 
♦  (2“  (-k) )  , {k , n , M- 1)3  . {n , 0  ,M- 1)3 
♦Exp [- (S*y/ (2*M) ) 3 / (2“M) , 

<y, 0.04283, 0.97772, 0.02)3 ; 
p2*Table [{y 1 [ [i3  3 ,pl [ [i3  3  *y 1 1 [ [i3  3 ) , 

{i, 1, Length [yl3)3 ; 

p3*Interpolation[p2,  InterpolationOrder  ->  13 ; 
p4-NIntegrate [p3 [y3 , {y , 0 . 04283 , 0 . 96283)3 ; 

Print [p433 


3.3  Convolutional  coding 


P«[ . 3;  . 

len-length(p) ; 

M-[4,3,2,3, 15, 104,540, 2520, 110483; 


r«M(2)/M(l);  . 

Pjl-C . 3;  . 

for  d«0:4  . 

pj-°; 

pj*M(4)+d; 

if  rem(j,2)**l  . 

for  i«(j+l)/2:j 

pj»pj+bino(j ,i)*p.‘i.*(l-p) .‘(j-i) ; 

end 


Input  p 

. Matrix  order  as 

n,k,v,df ,wd,wd+l,wd+2, . . . 

Solve  for  convolution  code  in 
matrix 

Solve  five  Pj  values 
d  i . e . ,  df  which  value  depend  on 
weight  wd+4 


Check  odd  or  even  case 
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else 

for  i«j/2+l:j 

p j  «p j+bino ( j , i ) *p . “ i . * ( l -p) . “ ( j -i ) ; 

end 

p j*pj+0 . 5*bino ( j , 0 . 5*j ) *p . * (0 . 5* j ) . * (1-p) . “ (0 . 5* j ) ; 

end 

pjl-[pjl;pj]  ; 

end 

pb»0; 

for  h«l:5 

pb*pb+M(h+4) *pj 1 (h, : ) ; 
end 

pb*pb/M(2) ; 

plot(s.pb)  . Plot  statement 
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